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Abstract 

The transport process in the fluid phase dur- 
ing the growth of a crystal has a profound influ- 
ence on the structure and quality of the solid 
phase. In vertical growth techniques the fluid 
phase is often subjected to vertical temperature 
and concentration gradients. The main objective 
of the present work is to obtain more experimental 
data on convective flows in enclosures with verti- 
cal temperature or concentration gradients. 

Among actual crystal systems the parameters 
vary widely. The parametric ranges studied herein 
for mass transfer experiments are mainly dictated 
by the electrochemical system employed herein to 
impose concentration gradients. Temperature or 
concentration differences are maintained between 
two horizontal end walls. The other walls are 
kept insulated. 

Experimental measurements and observations 
were made of the heat transfer or mass transfer, 
flow patterns, and the mean and fluctuating tem- 
perature distribution. For purely thermal cases, 
the experiment covers the ranges of thermal 
Rayleigh number Ra - 9x 1 0 6 to 2xlO y , aspect ratio 
Ar = 0.65 to 2.6, Prandtl number Pr = 4.0 to 7.0. 
The method used to visualize the flow pattern in 
the thermal cases is an electrochemical ph-indica- 
tor method. For purely solutal cases the experi- 
ment covers the ranges of solutal Rayleigh number 
Ra m = l.OxlO 9 to 2.4xlO n , Ar * 0.68 to 1.34, 

Sc = 2000 to 2100. Laser shadowgraphs is employed 
to visualize flow patterns In the solutal cases. 

The results of the present experiments for 
the purely thermal cases are adequately correlated 
by the equation Nu = 0.0862 Ra 0,322 in the range 
9x1 o 6 < Ra < 2xl0 9 . For the purely solutal cases, 
the mass transfer data can be correlated by the 
equation Sh * 0.1326 (Sc Gr m )°- 328 in the range 
IxlO 9 < Ra m < 2.4x10*’ . 

Nomenclature 

Ar aspect ratio (=d/L) 

Cfc bulk concentration of CUSO4 

C] constant of heat transport correlation 

C2 constant of mass transport correlation 

d height of enclosure 
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D solutal diffusivity 

F Faraday’s constant 

Gr thermal Grashof number 

Gr m solutal Grashof number 

g gravitational acceleration 

h average mass transfer coefficient 

I current 

ij_ limiting current density 

L length of enclosure' 

N buoyancy ratio 

Nu Nusselt number 

n exponential constant in transport correla- 

tion 

Pr Prandtl number 

Ra thermal Rayleigh number (=Gr Pr) 

Ra m solutal Rayleigh number (=Gr m Sc) 

Sc Schmidt number 

Sh Sherwood number 

T temperature 

Th hot wal 1 temperature 

Tc cold wall temperature 

t time 

V voltage 

H depth of enclosure 

x,y,z coordinates 

x,y,z dimensionless coordinates 

z* charge number of copper 

g volumetric coefficient of solutal expan- 

sion 

[3 volumetric coefficient of thermal expan- 

sion 

6 S thickness of solutal boundary layer 


V 


kinematic viscosity 
a thermal diffuslvity 

Introduction 

The growth of crystals is a coupled process 
of heat and mass transfer, fluid flow, phase trans- 
formations, and chemical reactions. Although the 
technology of crystal growth is well grounded in 
physical chemistry, the treatment of the transport 
phenomena has been relatively rudimentary and 
empirical. Recently, it was recognized that con- 
vection, as a result of temperature or concen- 
tration gradients , 1 occurs in many aspects of 
materials processing including solidification, oxi- 
dation of surface material, and crystal growth. 

The rate of crystal growth Is often limited by 
the mass transfer because it is the slowest proc- 
ess. In almost all the commonly used techniques 
the parent phase is a fluid. As in most nonequi- 
librium processes involving a fluid, convection 
plays a dominant role in that It affects the fluid- 
phase composition and temperature at the phase 
interface. The fluid motion in many cases Is 
determined by the locally generated motion. The 
flows are due either to the effects of temperature 
and/or concentration gradients in a body-force 
field . 2 In some crystal -growth techniques, the 
completely confined fluid phase is subject to ver- 
tical temperature and/or concentration gradients. 

To gain insights on such flows, an experimental 
program was initiated to study flows in rectangu- 
lar enclosures with vertical temperature or concen- 
tration gradients between horizontal end walls. 

The present work for purely thermal cases is 
intended to give a comprehensive understanding of 
natural convection heat transfer in enclosures with 
thermally insulated vertical walls by systemati- 
cally varying the parameters Ra and Ar. The 
fluid is tap water in these thermal experiments. 

An electrochemical system based on a diffusion- 
controlled electrode reaction was employed to cre- 
ate the vertical concentration gradients. The 
apparatus was a rectangular cavity with the hori- 
zontal end walls being the electrodes. All other 
walls were insulated. The fluid was a copper sul- 
phate (CuS 04 )-acid (H2SO4) solution for purely sol- 
utal experiments. 

Experimental Design 
Dimensionless Parameters 

Based on the basic differential equation for 
natural convection in an enclosure with thermal or 
solutal buoyancy forces (Fig. 1), it can be shown 
that the following dimensionless parameters are 
important in the present work. 

Thermal Grashof number, Gr = g 3 A Td^/v 2 
Prandtl number, Pr - v/a 
Schmidt number, Sc * v/D 
Solutal Grashof number, Gr m - g (5 A Cd^/v 2 
Aspect ratio, A * d/L 

With the present system the ranges of the 
dimensionless parameters covered for purely ther- 
mal cases are Pr » 4.0 to 7.0, Gr * 1.3xl0 6 to 
5.0x10°, Ar = 0.65 to 2.6, and for purely solutal 
cases are Sc - 2000 to 2100, Gr m • 5.0x10$ to 


1.2x10®, Ar * 0.68 to 1,34. All fluid properties 
are evaluated at the average temperature of the 
hot and cold walTC in the thermal cases. In the 
solutal cases fluid properties are evaluated at 
the average concentration of the bulk solution. 

Test Apparatus 

A sketch of the experimental system is pre- 
sented in Fig. 1 . The test cell is a rectangular 
enclosure formed with four Insulating pi ex i g 1 as 
plates and two copper plates (both thermally and 
electrically conductive). The width, L, is 7.7 cm 
or 8.2 cm and the height is variable so that a 
range of aspect ratios can be covered. The depth 
of the enclosure is relatively large (23.9 cm or 
21.2 cm). The two copper horizontal end walls are 
0.7 cm thick. 

For the purely thermal cases, an electrical 
heating mat is bonded to the back of the bottom 
copper plate and a circulating water system cools 
the top copper plate (Fig. 2(a)), so that a uniform 
temperature can be imposed on the top surface. The 
plexiglas surfaces are wrapped with 5 cm thick 
fiberglass insulation to ensure an adiabatic bound- 
ary condition. There are six thermocouples embed- 
ded In each copper plate to measure the uniformity 
of the temperature over the entire surface. A 
thermocouple probe with a two-dimensional travers- 
ing mechanism is inserted into the test section to 
measure the temperature distribution inside the 
enclosure. A Leeds and Northrup millivolt potenti- 
ometer is used .to measure the thermoelectric emf 
from the copper-constantan thermocouples. 

The test field is a copper-sulphate-acid solu- 
tion (CUSO4 + H 2 SO 4 + H 2 0 ) for the purely solutal 
cases. When a voltage is applied to the elec- 
trodes. copper dissolves into the solution at the 
anode and is deposited at the cathode (Fig. 2(b)). 
As a result, the density of the fluid near the 
cathode (anode) becomes lower (higher) than that 
of the bulk of the solution. The migration of 
cupric Ions in the electric field is eliminated by 
adding sulphuric acid to the solution, which acts 
as a supporting electrolyte, and thus the transport 
of the cupric ions is controlled only by diffusion 
and convection. The equipment for measuring the 
current and potential are shown schematically in 
Fig. 2. A dc power supply was used in the tests 
as the power source to the test cell and two digi- 
tal multimeters are necessary to obtain values of 
current correspond! ng to each potential setting. 

A potential recorder was used to record the volt- 
age as a function of time. 

Test Procedure 

The method used to visualize the flow patterns 
in the present work for purely thermal cases is an 
electrochemical method called the ph-indlcator 
method. ^ For this test thymol blue Is added to 
the water in the container to produce 0.01 percent 
by weight solution, and then the ph of the solu- 
tion Is adjusted to place it on the yellow side of 
the end-point. Two electrodes are placed in the 
solution. When a dc voltage (about 6 V in the 
present experiment) is applied between the elec- 
trodes, a small volume of colored solution forms 
around one of the electrodes and moves from the 
electrode with the fluid motion forming a neu- 
trally buoyant marker. The color change occurs 
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because the color of the basic form of thymol blue 
(blue) is different from that of the acid form 
(yel low) . 

Although the temperature of the copper walls 
can be easily measured by thermocouples, the con- 
centration levels at the walls can not be so eas- 
ily determined. One relatively simple way to set 
the concentration level at the cathode in the 
present system is to adjust the cell potential in 
such a way that the saturation (limiting) current 
is obtained. Under the limiting current condition 
the ion concentration at the cathode surface is 
zero. In other words, the change in concentration 
across the solutal boundary layer along the cath- 
ode ((AC) ca thode ) 1s c b- The limiting-current 
condition is valid at the cathode but not along 
the anode surface. However, since the net mass 
fluxes at the cathode and anode are considered to 
be equal and since the concentration outside the 
solutal boundary layers of both walls is Cb, it 
' is reasonable to expect that the average change in 
concentration across the solutal boundary layer 
along the anode <<AC) ar , 0 de ) 1s nedrl V ec l ua ] to 
(AC)cathode- Then, * he overall concentration dif- 
ference across the cell is AC = <AC) ca thode + 

(AC) a node * 2C b and thus Gr m is def1ned based 
on AC = 2Cb in the present work. The above con- 
cept assumes that Cb remains constant, but as 
convection develops in the cell the concentration 
level remains constant only on the average. More- 
over, due to electrolysis the copper wall surfaces 
become rough slowly with time, which causes the 
current densities at the walls to decrease and the 
cathode deviates gradually from the limiting- 
current condition. For this reason, the duration 
of each run for the pure solutal cases is limited 
to 2 hr, far shorter than the time required to 
attain steady solutal convection. It should thus 
be evident that the concentration boundary condi- 
tions are not as well-defined as the thermal 
conditions. After each run the electrodes are 
carefully cleaned and polished. 

The cell voltage is increased stepwise by a 
manual control and the corresponding total current 
is read after each steady state is attained. This 
procedure is repeated until a potential-current 
plateau is obtained. The potential stepping rate 
Is chosen as 20 mV/40 sec. The flow pattern and 
the concentration distributions are visually 
observed. Since the density of the solution 
changes with cupric ion concentration a qualita- 
tive picture of the flow configuration is obtained 
from the shadowgraph techniques. 

Experimental Results 
Purely Thermal Cases 

The flow structures were studied by the afore- 
mentioned ph-indlcator technique. Based on these 
observations the flow structures are shown in 
Figs. 3 to 5. "Thermals" rising from the bottom 
of the test cell and falling from the top of the 
test cell are evident. The dominant motion in the 
test cell is that of the thermals. 5 * 6 

Figure 6(a) and (b) shows typical mean temper- 
ature profiles. The mean temperature distributions 
are presented in the dimensionless form. The mag- 
nitude of the temperature, Tz - Tc, is normalized 
by the temperature difference, AT, across the 
layer. The subscript z denotes a condition at a 


distance z from the bottom plate. It was found 
that the temperature changed sharply near the top 
and bottom plates^'and the temperature distribution 
of the bulk flow was uniform. To record tempera- 
ture disturbances Induced by the thermals, the 
thermocouple probe was placed as near as possible 
to the bottom and top surfaces of the test cell. 
Time series records of the thermocouple output are 
shown in Figs. 7 and 8. According to Figs. 7 and 
8 the amplitudes of the temperature oscillations 
increase near the top and bottom plates. 

The Nusselt number data for water (Pr = 5.5 
to 7.0) for various values of Ar are presented 
in Fig. 9. The predictions of Nu by others b * b 
are also shown in this figure for comparison. Nu 
was found to be proportional to Ra 0 - 322 as seen 
in Fig. 9. The value of the exponent (n = 0.322) 
points to a heat-transfer mechanism by turbulent 
eddies. However, there appears to be a small but 
noticeable difference between the present data 
(Hr = L/W = 3.22) and others' data (Hr == 1). Nu 
seems to increase slightly with decreasing Hr 
(L/W). This dependence of Nu upon Hr may be 
explained if one considers the probably weak hori- 
zontal flow motion which can be detected in 
Figs. 3 and 4. 

Purely Solutal Cases 

In the present electrochemical system opera- 
ted under isothermal conditions the deposition of 
cupric ions on the cathode wall leaves behind less 
dense fluid. The lighter fluid rises toward the 
top of the test section. The downward flow in the 
boundary layer near the anode feeds heavier fluid 
toward the bottom wall. 7 Simple shadowgraph stud- 
ies (Figs. 10 and 11) suggested a very complex 
flow structure over the Sc Gr m range 1x10 to 
2 . 4x 1 0 1 1 . 

It may be useful to compare in detail the 
purely solutal convection flow structure with that 
for purely thermal convection because the overall 
flow structures look somewhat similar for both 
cases. An important fundamental contrast, how- 
ever, in the flow structures produced by purely 
solutal convection and purely thermal convection 
is illustrated in Figs. 11 and 5. In the heat 
transfer situation there are far fewer convection 
columns (thermals) than for mass transfer. 

The most convenient method for studying rates 
of mass-transfer by the free convection mechanism 
Is by measuring the maximum rate of the electrode 
reaction (i.e., limiting current) which may be 
maintained in the quasi-steady state. With the 
limiting current per unit area of the cathodes, 

1 Lf known the mass transfer according to Ref. 7 
can be calculated as 


h = 


C 


b 


The dimensionless mass transfer rate, the Sherwood 
number. Is defined as 


Sh = 


h d 

h D 


In the present electrochemical system the dif- 
fusion coefficient is very small so a long time is 
required to obtain truly solutal convection. Thus, 
the mass transfer rate measured herein In not a 
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steady-state value. However, the solutal boundary 
layers should be established on a time scale of 
S s / D , which Is on the order of 10 sec and Is small 
compared to the time at which the limiting current 
is measured. The measured mass-transfer rate Is 
not expected to differ much from the steady-state 
value. This accuracy is probably acceptable par- 
ticularly in view of the uncertainty in the concen- 
tration boundary conditions. The experimental 
error in Sh is estimated to be ±10 percent. A 
correlation for Sh, as a function of solutal 
Rayleigh number (Sc Gr m ) , is presented in Fig. 12. 

Summary of Resul ts 

The results of the present experiments indi- 
cate that for the range of aspect ratios and 
Grashof (or Rayleigh) numbers covered: 

1. The results for the purely thermal cases 
are adequately correlated by the equation 

Nu = 0.862 Ra 0 - 322 
in the range 9xl0 6 < Ra < 2xl0 9 . 

2. Temperature oscillations and the exponent 

0.322 found in the purely thermal cases reflect a 
high degree of flow instability in the present 
work . 

3. For the purely solutal cases the mass 
transfer data can be correlated by the equation 

Sh = 0.1326 (Sc Gr m ) 0 - 328 
in the range 1x10^ < Ra m < 2.4x10''. 

In the present work, comparison of the new 
correlations of mass-transfer data to the results 
of heat-transfer investigations could serve to 
give us some enlightenment about the mechanism of 
the transport process under discussion. 

The figures 9 and 12 show clearly that all 
results are well correlated by the equations: 

Nu = C] Ra n , 

Sh = C 2 Ra m ", 

from a practical point of view n - 1/3 Is accepta- 
ble as a regression parameter, at least within the 
experimental Ra (or Ra m ) range. 

It enables a general understanding of convec- 
tive flow for all important vertical crystal 
growth methods . 
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FIGURE 2. - EXPERIMENTAL SYSTEMS. 
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<C> (d) 

FIGURE 3. - FLOW PATTERNS NEAR THE HEATED LOWER PLATE (PURELY THERMAL CASE) WITH 
Ar = 0.65; Ra « I.SIxlO 7 ; t-j < = 5 SEC; 7 = 0 TO 0.5. 
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(c) t 3 . (d) t^. 


FIGURE M. - FLOW PATTERNS NEAR THE LOWER PLATE (PURELY THERMAL CASE) WITH Ar = 0.65; 
Ra = 1.51x10 7 ; t-j < t 2 < t 3 < tj,; At = 5 sec; y = 0.6 TO 1.0. 
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(a) t 1 



(b) t 2 = ti + 5 SEC. 

FIGURE 5. - FLOW PATTERNS NEAR THE COLD 
TOP PLATE (PURELY THERMAL CASE) WITH 
Ar = 1.30; Ra = 1.35x10 8 ; t 1 < t 2 ; 

At = 5 SEC. 
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FIGURE 8. - INSTANTANEOUS TEMPERATURE PROFILES AT x = 0.52. y = 0.17 FOR PURELY THERMAL CASE WITH Ar = 1.30, Ra = 1.35x10' 
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FIGURE 9. - CORRELATION OF Nu FOR PURELY THERMAL CASE, 
Hr = 0.322 (PRESENT CASE). 
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(b) t2 = *1 + 5 sec. 

FIGURE 10. - FLOW PATTERNS NEAR THE CATH- 
ODE LOWER PLATE (PURELY SOLUTAL CASE), 
Ar = 1.02, Ra m = 2.08X10 11 . ^ < t 2 * 

At = 5 SEC. 
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(b) t 2 = t-, + 15 sec. 


FIGURE 11. - FLOW PATTERNS NEAR THE ANODE 
TOP PLATE (PURELY SOLUTAL CASE). 

Ar = 1.02. Ra m = 2.08X10 11 . t, < t 2 * 

At = 15 SEC. 


10x10 2 



FIGURE 12. - CORRELATION OF Sh FOR PURELY SOLUTAL CASE. 
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